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Abstract. This work presents the results from a systematic search for evidence of temporal changes (i.e., non-stationarity) 
associated with spectral variations in 3C 390.3, using data from a two-year intensive RXTE monitoring campaign of this 
broad-line radio galaxy. In order to exploit the potential information contained in a time series more efficiently, we adopt a 
multi-technique approach, making use of linear and non-linear techniques. All the methods show suggestive evidences for non- 
stationarity in the temporal properties of 3C 390.3 between 1999 and 2000, in the sense that the characteristic time-scale of 
variability decreases as the energy spectrum of the source softens. Flowever, only the non-linear, “scaling index method” is able 
to show conclusively that the temporal characteristics of the source do vary, although the physical interpretation of this result is 
not clear at the moment. Our results indicate that the variability properties of 3C 390.3 may vary with time, in the same way as 
they do in Galactic black holes in the hard state, strengthening the analogy between the X-ray variability properties of the two 
types of object. This is the first time that such a behavior is detected in an AGN X-ray light curve. Further work is needed in 
order to investigate whether this is a common behavior in AGN, just like in the Galactic binaries, or not. 

Key words. Galaxies: active - Galaxies: nuclei - X-rays: galaxies 


1. Introduction 

Among the basic properties characterizing an Active Galactic 
Nucleus (AGN) the flux variability, especially at X-rays, is 
one of the most common. Previous multi-wavelength variabil¬ 
ity studies (e.g., Edelson et al. 1996; Nandra et al. 1998) have 
shown that AGN are variable in every observable wave band, 
but the X-ray flux exhibits variability of larger amplitude and 
on time-scales shorter than any other energy band, indicating 
that the X-ray emission originates in the innermost regions of 
the central engine. 

However, even though the discovery of X-ray variability 
dates back more than two decades, its origin is still poorly 
understood. Nevertheless, it is widely acknowledged that the 
variability properties of AGN are an important means of prob¬ 
ing the physical conditions in their emission regions. The rea¬ 
son is that using energy spectra alone it is often impossible to 
discriminate between competing physical models and thus the 
complementary information obtained from the temporal analy¬ 
sis is crucial to break the spectral degeneracy. 

Over the years, several variability models have been pro¬ 
posed, involving one or a combination of the fundamental com¬ 
ponents of an AGN: accretion disk, corona, and relativistic jet. 
These models can be divided into two main categories: 1) in- 
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trinsically linear models, as the shot noise model (e.g., Terrel 
1972), where the light curve is the result of the superposition of 
similar shots or flares produced by many independent active re¬ 
gions, and 2) non-linear models, which require some coupling 
between emitting regions triggering avalanche effects, as the 
self-organized criticality model (e.g., Mineshige et al.l994), 
or models assuming that the variability is caused by varia¬ 
tions in the accretion rate propagating inwards (e.g. Lyubarskii 
1997). The “rms-flux relation” recently discovered by Uttley & 
McHardy (2001), along with detection of non-linearity in dif¬ 
ferent AGN light curves (i.e., in 3C 345 by Vio et al. 1992; in 
3C 390.3 by Leighly & O’Brien 1997; in NGC 4051 by Green 
et al. 1999; in Ark 564 by Gliozzi et al. 2002), favors non-linear 
variability models, but there is no general consensus on the na¬ 
ture of the X-ray variability yet. 

Because of their brightness, the temporal and spectral prop¬ 
erties of Galactic black holes (GBHs) are much better known 
and can be used to infer information on their more power¬ 
ful, extragalactic analogs, the AGN. It is well established that 
GBHs undergo state transitions (see McClintock & Remillard 
2003 for a recent review), switching between “low”, “inter¬ 
mediate”, “high”, and “very high” states, which are all unam¬ 
biguously characterized by a specific combination of energy 
spectrum and power spectral density (PSD). Importantly, even 
when in the same state, the timing properties of GBHs may 
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vary with time, often accompanying changes of the spectral 
properties. For example, Pottschmidt et al. (2003) showed that 
changes in the PSD of Cyg X-1 are associated with changes 
in the slope of the energy spectrum. Specifically, during the 
“low/hard state”, the PSD characteristic frequencies increase, 
as the spectral slope becomes steeper. 

Transitions between “different” states have not been ob¬ 
served yet in AGN. The reason may be associated with the fact 
that their flux variability is characterized by much longer time- 
scales. Nevertheless, a correspondence between GBH spec¬ 
tral states and some AGN classes has been hypothesized. 
Specifically, Broad-Line Seyfert 1 galaxies have been associ¬ 
ated with the GBH low/hard state, and Narrow-Line Seyfert 1 
galaxies with the GBH high/soft state, (for a recent discussion 
on this topic see, for example, McHardy et al. 2004). 

Although genuine spectral transitions in AGN are unlikely 
to be detected due to the longer time-scales involved, it is per¬ 
haps possible to detect changes in the temporal properties (i.e., 
loss of stationarity) associated with spectral changes, similarly 
to what observed in GBHs within the same spectral state. To 
this end, of crucial importance are the numerous monitoring 
campaigns carried out by RXTE in the last decade, which have 
produced data sets suitable for searching for non-stationarity in 
the AGN light curves, where long-term spectral changes proved 
to be ubiquitous. 

In this work, we present the results from a systematic 
search for evidence of temporal changes associated with sig¬ 
nificant spectral variations in the case of the broad-line radio 
galaxy (BLRG) 3C 390.3. We use data from a two-year in¬ 
tensive RXTE monitoring campaign. Previous X-rays studies 
with ROSAT, and RXTE (Leighly & O’Brien 1997, Gliozzi et 
al. 2003) demonstrated that 3C 390.3 is one of the most vari¬ 
able AGN on time-scales of days and months. Our motivation 
stems from the results of a long monitoring campaign carried 
out with RG^ATHRI, which showed evidence for possible non- 
stationarity (and non-linearity) in the soft X-ray light curve of 
3C 390.3 (Leighly & O’Brien 1997). 

2. Observations and Data Reduction 

We use archival RXTE data (PI; Leighly) of 3C 390.3. This 
source was observed by RXTE for two consecutive monitoring 
campaigns between 1999 and 2001. The first set of observa¬ 
tions was carried out from 1999 January 8 to 2000 February 
29, and the second one from 2000 March 3 to 2001 February 
23. Both campaigns were performed with similar sampling: 
3C 390.3 was regularly observed for ~ 1000-2000 s once ev¬ 
ery three days. The observations were carried out with the 
Proportional Counter Array (PCA; Jahoda et al. 1996), and the 
High-Energy X-Ray Timing Experiment (HEXTE; Rotschild 
et al. 1998) on RXTE. Here we will consider only PCA data, 
because the signal-to-noise of the HEXTE data is too low for a 
meaningful timing analysis. 

The PCA data were screened according to the follow¬ 
ing acceptance criteria: the satellite was out of the South 
Atlantic Anomaly (SAA) for at least 30 minutes, the Earth 
elevation angle was > 10°, the offset from the nominal op¬ 
tical position was < 0?02, and the parameter ELECTRON-2 


was <0.1. The last criterion removes data with high particle 
background rates in the Proportional Counter Units (PCUs). 
The PCA background spectra and light curves were deter¬ 
mined using the L7 - 240 model developed at the RXTE 
Guest Observer Eacility (GOP) and implemented by the pro¬ 
gram pcabackest v.2.1b. Since the two monitoring cam¬ 
paigns span three different gain epochs, the appropriate files 
(pcaJ}kgd_cmfaintn_e3(4,5)bv2QQ3ll23.mdl, depending on 
the observation date), provided by the RXTE Guest Observer 
Eacility (GOP), were used to calculate the background light 
curves. This model is appropriate for “faint” sources, i.e., those 
producing count rates less than 40 s ' PCU“'. All the above 
tasks were carried out using the FTOOLS v.5.1 software pack¬ 
age and with the help of the rex script provided by the RXTE 
GOP, which also produces response matrices and effective 
area curves appropriate for the time of the observation. Data 
were initially extracted with 16 s time resolution and subse¬ 
quently re-binned at different bin widths depending on the ap¬ 
plication. The current temporal analysis is restricted to PCA, 
STANDARD-2 mode, 2-20 keV, Layer 1 data, because that 
is where the PCA is best calibrated and most sensitive. Since 
PCUs 1,3, and 4 were frequently turned off (the average num¬ 
ber of PCUs is 2.25 and 1.95 in the 1999 and 2000 observa¬ 
tions, respectively), only data from the other two PCUs (0 and 
2) were used. All quoted count rates are for two PCUs. 

The spectral analysis of PCA data was performed using the 
XSPEC V. 11 software package (Arnaud 1996). We used PCA 
response matrices and effective area curves created specifically 
for the individual observations by the program pcarsp, taking 
into account the evolution of the detector properties. All the 
spectra were rebinned so that each bin contained enough counts 
for the;^^^ statistic to be valid. Pits were performed in the energy 
range 4-20 keV, where the signal-to-noise ratio is the highest. 

3. The X-ray Light Curve 

Pigure[2shows the background-subtracted 2-20 keV count rate 
and hardness ratio (HR=7-20 keV/2-5 keV) light curves of 
3C 390.3. Pilled circles represent data points from the first 
monitoring campaign, whereas open diamonds are used for the 
second campaign. A visual inspection of the HR and count 
rate light curves indicates the presence of a typical Seyfert- 
like spectral evolution (the energy spectrum softens when the 
source brightens), and suggests the existence of different vari¬ 
ability patterns during the two monitoring campaigns. In the 
first one, 3C 390.3 exhibits an initial smooth decrease in the 
count rate (a factor ~2 in 100 days) followed by an almost 
steady increase (a factor ~3.5 in 100 days). After day ~ 200 
(from the beginning of the observation), the source shows 
significant, low-amplitude, fast variations. During the second 
monitoring campaign, the source does not show any long-term, 
smooth change in count rate, but only a flickering behavior 
with two prominent flares. This qualitative difference is asso¬ 
ciated with marked differences in the two mean count rates, 
6.4 + 0.2 s ' and 9.0 + 0.1 s ', and in the fractional variabil¬ 
ity amplitudes, Evar.gg = (33 ± 3)% and Evar.oo = (19 ± 2)%, 
calculated over each campaign (see, e.g., Vaughan et al. 2003 
for a detailed description of Eyar and the estimate of its uncer- 
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Fig. 1. X-ray light curves of the background-subtracted count rate in the 2-20 keV band (top panel) and of the X-ray color 7-20 
keV/2-5 keV (bottom panel) from RXTE PCA observations of 3C 390.3. Time bins are 5760s (~ 1 RXTE orbit). Filled circles 
represent data points from the hrst monitoring campaign from 1999 January 8 to 2000 February 29, whereas open diamonds are 
data points from the second monitoring campaign from 2000 March 3 to 2001 February 23. The dashed lines in the top panel 
indicate the intervals used for the time-resolved spectral analysis (see §9). 


tainty). The uncertainties have been computed assuming that 
the data distribution is normal and the data are a collection of 
independent measurements. However, the data points of AGN 
light curves are not independent but correlated and thus the ap¬ 
parent difference between the mean count rates and Fy^ could 
simply result from the red-noise nature of the variability pro¬ 
cess, combined with its randomness (see, e.g., Vaughan et al. 
2003). 

4. Stationarity Tests 

Almost all methods employed in timing analysis (e.g., power 
spectrum, structure function, auto-correlation analysis) assume 
some kind of stationarity. Generally speaking, a system is con¬ 
sidered stationary, if the average statistical properties of the 
time series are constant with time. More specihcally, the sta¬ 
tionarity is dehned as “strong” if all the moments of the proba¬ 
bility distribution function are constant, and “weak” if only the 
moments up to the second order (i.e., mean and variance) are 
constant. Usually, non-stationarity is considered an undesired 
complication. However, in some cases the non-stationarity pro¬ 
vides the most interesting information, as in the case of GBHs 
switching from a spectral state to another. 

In previous temporal studies of AGN, no strong evidence 
for non-stationarity has been reported. Indeed, AGN light 
curves are considered to be “second-order” stationary, due to 
their red-noise variability. Here, we adopt several complemen¬ 


tary methods to investigate the issue of the stationarity in the 
3C 390.3 light curve. First, the temporal properties are stud¬ 
ied with traditional linear techniques in the Fourier domain, 
by estimating the PSD, and in the time domain with the struc¬ 
ture function analysis. Secondly, we probe directly the proba¬ 
bility density function associated with the time series. Finally, 
we utilize the scaling index method, a technique borrowed from 
non-linear dynamics, useful for detecting weak signals in noisy 
data. 

All the above methods are separately applied to the 1999 
and 2000 monitoring campaigns. We quantify the differences 
between the estimated functions (like the power spectrum, 
the structure function, etc.) using appropriate “statistics”. The 
probability that these statistics will have the observed val¬ 
ues are rigorously estimated with extensive Monte Carlo sim¬ 
ulations. These simulations allow us to assess whether the 
observed values are due to measurement uncertainties (i.e., 
Poisson noise), to the light curve sampling, to the red-noise 
nature of the variability process, or to intrinsic variations of the 
statistical properties of the source. 

4.1. Analysis in the Fourier Domain: Power Spectral 
Density 

We followed the method of Papadakis & Fawrence (1993) to 
compute the power spectrum of the 1999 and 2000 light curves, 
after normalizing them to their mean. We used a 3-day bin- 
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ning (i.e., equal to the typical sampling rate of the source) for 
both light curves. As a result, the light curves that we used are 
evenly sampled, but there are a few missing points in them; 30 
out of 140 and 15 out of 120 for the 1999 and 2000 light curves, 
respectively. These points are randomly distributed over the 
whole light curve, and we accounted for them using linear in¬ 
terpolation between the two bins adjacent to the gaps, adding 
appropriate Poisson noise. 

Figure |2] shows the PSD of the 1999 and 2000 light curves 
(filled circles, and open diamonds, respectively). Each point in 
the estimated PSD represents the average of 10 logarithmic 
periodogram estimates. We fitted both spectra with a simple 
power law (P(/) k / “). The model provides a good descrip¬ 
tion of both PSDs; the best fitting slope values (shown in Fig|2j 
are consistent within the errors with the value of 2 + 0.3 in both 
cases. 

A closer look at the PSD plot indicates a possible low fre¬ 
quency flattening of the 2000 PSD between log(v) = -6.5 and 
-7, corresponding to a time interval between ~35 and 115 days. 
In order to investigate this issue further we split the two lowest 
frequency points into two (to increase the frequency resolution) 
and refitted the 2000 PSD with a broken power-law model. We 
kept the slope below and above the frequency break fixed to 
the values 1 and 2, respectively, and kept as free parameters 
the break frequency and normalization of the PSD. The model 
provides a good fit to the 2000 PSD, but statistically not better 
than a simple power law according to an F-test. The best fitting 
break frequency is 2.5 x 10“^ Hz, which corresponds to a break 
time-scale of ~ 46 days. The 90%-confidence lower limit to 
this time-scale is 33 days, while, due to the limited coverage of 
the observed PSD at lower frequencies, the upper limit is un¬ 
defined. We also tried to fit the 1999 PSD with this model, but 
without success: the best fitting break frequency turned out to 
be lower than the lowest frequency sampled by the 1999 light 
curve, indicating that there is no indication of any breaks in the 

1999 PSD. This power spectrum follows a power-law shape 
with slope a — 2 down to the lowest sampled frequencies. 

In order to quantify the comparison between the 1999 and 

2000 PSDs in a statistical way, we used the “5” statistic as de¬ 
fined by Papadakis & Lawrence (1995). To avoid the constant 
Poisson noise level, we considered the periodogram estimates 
up to frequencies ~ 7 x 10“^ Hz. We find that S = -1.3, which 
implies that the two light curves are not far away from the hy¬ 
pothesis of stationarity. Therefore, the differences suggested in 
the previous paragraph by the broken power-law model fitting 
results are not supported by this more rigorous statistical test. 
Light curves with denser sampling, which would result in a 
much higher frequency resolution at lower frequencies (where 
the potential differences between the two PSDs are observed), 
would be necessary for a more sensitive comparison between 
the 1999 and 2000 power spectra. 

Recent results from the power spectrum analysis of com¬ 
bined XMM-Newton and RXTE light curves of radio-quiet AGN 
(Markowitz et al., 2003; McHardy et al. 2004) show that the 
2-10 keV power spectrum of AGN follows a power-law 
form with a slope of 2 down to a “break” frequency, below 
which it flattens to a slope of 1. This break time-scale appears 
to scale with the black hole mass according to the relation: 



Fig. 2. RXTE PCA power density spectra using based on the 
1999 and 2000 light curves with 3-day bins, normalized to 
their mean (filled circles and open diamonds, respectively). The 
solid and dashed lines show the best fitting power law models, 
including the effects of the experimental Poisson noise level. 


fir(days) = Mbh/10®'^ Mq (Markowitz et al. 2003). Using the 
reverberation mapping technique, Peterson et al. (2004) esti¬ 
mate the black hole mass of 3C 390.3 to be Mbh ~ 2.9 x 10* 
Mq. In this case, the characteristic break time-scale of 3C 390.3 
should be ~ 90 days, which corresponds to a break frequency 
of ~ 10“^ Hz. This break frequency is entirely consistent with 
the broken power-law model fitting results in the case of the 
2000 light curve. Unfortunately, though, even the combined 

1999 and 2000 PSD does not extend to frequencies low enough 
for this break frequency to be firmly detected. 

It is worth noticing that the value assumed for the mass of 
the black hole harbored by 3C 390.3 is still a matter of debate. 
The value used in this work is the most recent and it is consis¬ 
tent within the errors with the mass estimate given by Kaspi et 
al. (2000). However, based on another reverberation mapping 
experiment on this source, Sergeev et al. (2002) find a longer 
reverberation lag and determine a higher mass of 2 x 10® Mq 
using a different method. In that case, the putative break fre¬ 
quency in the PSD would be located at a frequency of ~ 1/yr, 
which could not be probed with the current data set. 

We conclude that the combined 1999H-2000, 2 - 20 keV 
PSD of 3C 390.3 is well fitted by a simple power-law model 
at all sampled time-scales from 3“' - 100“* days“*. The same 
model fits well the individual 1999 and 2000 PSDs as well. 
There is an indication of a low-frequency PSD break (at a fre¬ 
quency which corresponds to a time-scale of ~50 days) in the 

2000 PSD, which is absent in the 1999 PSD. However, due to 
the low resolution of the power spectra at hand, we cannot be 
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Fig. 3. Structure function of 3C 390.3 during the 1999 (filled 
circles) and 2000 (open diamonds) monitoring campaign. Time 
bins are 3 days. For sake of clarity, error-bars were plotted only 
for the 1999 SF. 

certain of the significance of the difference between the two 
PSDs. 

4.2. Analysis in the Time Domain: Structure Function 

Next, we performed an analysis based on the structure func¬ 
tion (e.g., Simonetti 1985), a linear method which works in the 
time domain and has the ability to discern the range of time- 
scales that contribute to the variations in the data set. In prin¬ 
ciple, the SF should provide the same information as the PSD, 
as both functions are related to the auto-covariance function of 
the process R{t). The power spectrum is defined as the Fourier 
transform of /?(t), while the intrinsic SF at lag r is equal to 
2cr^ - R(t), where cr^ is the variance of the process (Simonetti 
et al., 1985). In fact, if the PSD follows a power law of the 
form P{f) oc (where / is the frequency), then SF oc 
(Bregman et al., 1990). 

We computed the structure functions for the first and sec¬ 
ond monitoring campaigns. The results are plotted in Fig|3 the 
filled circles refer to the 1999 light curve, whereas the open 
diamonds describe the 2000 SF. Both SFs have a power-law 
shape of the form SF oc t* , which is consistent with the results 
of the PSD analysis, since P(f) oc (see §4.1). Furthermore, 
the 2000 SF shows a flattening to SF oc t“, at a characteristic 
time lag of ~ 4 - 6 x 10^ s, followed by an oscillating be¬ 
havior. This value corresponds to a characteristic time-scale of 
~ 50-70 days, which is fully consistent with the putative break 
suggested by the 2000 PSD and corresponds approximately to 
the value expected on the basis of the black hole mass estimate 
(see discussion on §3.1). 

The 1999 and 2000 SFs look qualitatively different. 
However, the statistical significance of this difference cannot 
be assessed directly, since the the SF points are correlated and 



Fig. 4. Comparison between Dsp.obs (thick continuous line) and 
the distribution of Dsp.synth obtained with red-noise simula¬ 
tions. 

their true uncertainty is unknown. The errors shown in Fig. |3 
are representative of the typical spread of the points around 
their mean in each time-lag bin, hence they depend mainly on 
the the number of points that contribute to the SFs estimation at 
each time lag. To investigate in a quantitative way whether the 
difference between the two SFs can be simply ascribed to the 
intrinsic, red-noise randomness of the light curves or it actu¬ 
ally reveals a non-stationary behavior, the use of Monte Carlo 
simulations is necessary. 

For this reason, we performed a model-dependent numer¬ 
ical experiment assuming that both the 1999 and 2000 light 
curves are realizations of a process with mean count rate equal 
to the average value of the 1999 and 2000 light curves, and 
with a power spectrum of power-law index 2 down to a break 
frequency below which the power-law index changes to 1. We 
performed the experiment considering various PSD break time 
scales: fpr = 50,100,200, and 300 d. Here we present the only 
the results from fpr =300 d because they provide the most 
conservative limits on the significance of the observed differ¬ 
ences. In all cases, the PSD normalization was fixed in such 
a way that the model PSD will have the same value as the ob¬ 
served 1999H-2000best fitting power law at 10“^ Hz. Assuming 
this power spectrum model, we created two sets of 10^ syn¬ 
thetic light curves (using the method of Timmer & Koenig 
1995). Importantly, for each synthetic light curve, the appropri¬ 
ate RXTE background light curve (as provided by the FTOOL 
PCABACKEST) was added to each light curve, and then the ap¬ 
propriate Poisson noise was added to each point in them. The 
length of the synthetic light curves was 20 times longer than the 
length of the observed light curves, and their points were sep¬ 
arated by intervals of 160 s. In this way, we could simulate the 
effects of the red noise leakage in the estimation of the SF, al¬ 
though the aliasing effects were not accounted for. This should 
not be a serious problem though, as the steep power spectral 
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slope (-2) implies that these effects should not be very strong. 
Each synthetic light curve was then re-sampled in such a way 
to have the same length and sampling pattern as the observed 
light curves. 

It is worth noticing that, after fixing the PSD normaliza¬ 
tion, the synthetic light curves are not renormalized to the ob¬ 
served variance, but their variance is free to vary, allowing a 
larger scatter in the low-frequency power and hence maximiz¬ 
ing the effects of long-term flux variations. For completeness, 
we also performed a second set of “standard” red-noise simula¬ 
tions with synthetic light curves renormalized in such a way to 
have the same sampling pattern, mean, and variance as the 1999 
and 2000 light curve, respectively. No significant differences in 
the red-noise tests of SF are found using the two different nu¬ 
merical approaches. 

For each pair of synthetic light curves, we computed their 
SFs, and used the sum of the squared ratios between them at 
each lag, i.e. Dsf - 'ZISF 2 (Ti)/SFiin)]^, to quantify their 
difference. Based on the distribution of the Dsp.synth values (see 
Fig.|4j, we find that the probability that the Dsp.obs value would 
be as large as observed, is of the order of 44%, Therefore, the 
apparent differences between the observed 1999 and 2000 SFs 
shown in Fig.|3can be the result of the red-noise character of 
the variability process. 

However, we noticed that the synthetic SFs show the largest 
differences at lags larger than ~ 100 days (i.e., at those lags 
where the wiggles and bumps in the observed SFqo appear). 
In fact, when we restricted the computation of the Dsf up 
to lag=90 days (thus including the flattening of the observed 
S Fqq, but excluding the part which shows the highly erratic 
behavior at large lags) we found that probability of the differ¬ 
ence between the observed 1999 and 2000 SFs being due to 
red-noise effect decreases to 7.5%. 

The SF analysis results are similar to those from the power 
spectrum analysis. The SFqq shows a flattening around lag 
~ 50 - 70 days, which is absent in the S Fgg, suggesting that 
the power spectrum break may have indeed moved at higher 
frequencies during the 2000 observations, in agreement with 
the 1999 and 2000 PSDs. However, just like with the observed 
PSDs, we cannot be certain of the reality of this effect. 

4.3. Time-independent Analysis: Probability Density 
Function 

After utilizing methods which are based on the temporal or¬ 
der of the data points, here we investigate the stationarity of 
3C 390.3 with a method independent of the temporal order of 
the data points, the probability density function (PDF), which 
is the histogram of the different count rates recorded during 
the monitoring campaign. Fig.|3shows the normalized distribu¬ 
tions (r-(r))/(r) of the 1999 (thick solid line) and 2000 (dotted 
line) monitoring campaigns. The first transformation (r - (r)) 
brings both light curves to zero mean, the second, normalizes 
the amplitudes to unit mean. In this way, the relative amplitudes 
of the light curves are directly probed, and the PDF results can 
be directly compared to the power spectral ones, which use ex¬ 
actly the same convention. 



-1.0 -0.5 0.0 0.5 1.0 

(r-<r>)/<r> 


Fig. 5. Probability density function of count rates scaled and 
normalized with respect to the average count rate (r), during 
the 1999 (thick solid line) and 2000 (dotted line) monitoring 
campaign. The time bins are 32 s wide. 


For a better visual comparison, we have used small time 
bins of 32 s in Fig. |3 However, for a quantitative compari¬ 
son between PDFs (see below), we have used larger time bins 
to prevent any contamination from the Poisson noise. There 
is a clear qualitative difference between the two PDFs, with 
the 1999 distribution characterized by two broad peaks and the 
2000 PDF having a bell-like shape. This qualitative difference 
is confirmed by a Kolmogorov-Smirnov (K-S) test, which gives 
a probability of 99.99% that the two distributions are different. 

One may ask what are the effects on the PDF properties of 
the light curve bin size, whether the PDF is really time inde¬ 
pendent (a question that is directly related to the length of the 
light curve), and whether the K-S test is appropriate to quantify 
the statistical difference between two PDFs or not. The answers 
to these questions are not clear. In all cases, the most important 
issue is the red-noise character of the AGN light curves (in all 
wavebands). As the size of the light curve bins decreases, the 
number of the points in the PDF increases, and the sampled 
PDF is better defined. However, as the points added are heav¬ 
ily correlated due to the red-noise character of the variations, 
no “extra” information is added to the PDF. On the other hand, 
since the “sensitivity” of the K-S test increases with the number 
of points in the sampled PDF (which are assumed to be inde¬ 
pendent), there is a possibility that differences which are purely 
due to the stochastic nature of a possible stationary process will 
be magnified and the K-S test will erroneously indicate “signif¬ 
icant” differences between two PDFs. Finally, one should also 
consider the effects of the Poisson noise to the sampled PDFs. 
In the cases where the source signal has a much larger ampli¬ 
tude than the variations caused by the experimental noise (as is 
the case here), we should not expect this effect to be of any sig¬ 
nificance. In any case, the Poisson noise should produce similar 
PDFs (once they are normalized to their mean) so any signifi¬ 
cant differences cannot be due to this effect. 
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Fig. 6. Comparison between DpDp.obs (thick continuous line) 
and the distribution of OpDF.synth obtained with red-noise simu¬ 
lations. 

To assess the apparent PDFs differences shown in Fig. |3 
in a more rigorous way, we have used the 10 ^ synthetic light 
curves with tbin=160 s, obtained from the Monte Carlo simu¬ 
lations described in §4.2. For each pair of synthetic PDFs, we 
computed the sum Dpop - YlPDF\{i) - PDF 2 (iy\^, in order 
to quantify their differences. A comparison of the distribution 
of the 10 ^ synthetic DpDp.synth with the observed value DpDp.obs 
(represented by a thick solid line in Fig.|^ reveals that a value 
as large as Dpop.obs appears in 14.3% of all the synthetic light 
curves. Consequently, the apparent difference between PDF 99 
and PDFoo could be caused by red-noise effects. 

Looking closer at the synthetic PDFs that show the largest 
Dpop values, we noticed that they are caused mainly by differ¬ 
ences at \{r-{r))!{r)\ > 0 / 6 , where there is only a small fraction 
of the total number of points in the histograms. For that reason, 
we divided the x-axis into three sub-intervals (- 1 /- 0 . 6 ; - 0 . 6 / 0 . 6 ; 
and 0.6/1) and calculated Tipop in each one of them. When we 
consider the middle interval, the probability that the difference 
between the observed PDFs can be explained to red-noise ef¬ 
fects is reduced to 12 . 6 %. 

In summary, the 1999 PDF appears to be broader than the 
2000 PDF, which shows a well defined, uniform distribution 
around its mean. It shows a bimodal shape, with a large frac¬ 
tion of points having a count rate smaller than average. This is 
qualitatively consistent with the results from the previous sec¬ 
tions. For example, in the case of a PSD break at low frequen¬ 
cies, the observed variations are dominated by the components 
with long periods. If the PDF is estimated with data sets shorter 
than the PSD break time-scale, we would expect the points to 
show a broad distribution around the light curve mean, as in 
the case of the 1999 PDF. Indeed, our results in the case of a 
PSD break at 1/50 days ' show that the observed 1999 PDF is 
entirely inconsistent with this assumption, contrary to the ob¬ 
served 2000 PDF which is consistent. However, as before with 


the results from the PSD and SF analysis, a red-noise process 
with a PSD break at 1/300 days * is in rough agreement with 
both the 1999 and 2000 PDFs. So, we conclude that observed 
1999 and 2000 PDF differences possibly support the hypothe¬ 
sis of non-stationarity, but not at a high significance level. 

4.4. Non-linear Analysis: Scaling Index Method 

Non-linear methods are rarely employed in the analysis of 
AGN light curves, partly because these methods are less de¬ 
veloped than the linear ones, partly because it is implic¬ 
itly assumed that AGN light curves are linear and stochastic. 
However, non-linear methods can be useful not only to char¬ 
acterize the nature of chaotic deterministic systems, but also to 
discriminate between two time series, regardless of the fact that 
they are linear or non-linear. 

When methods of non-linear dynamics are applied to time 
series analysis, the concept of phase space reconstruction rep¬ 
resents the basis for most of the analysis tools. In fact, for 
studying chaotic deterministic systems it is important to es¬ 
tablish a vector space (the phase space) such that specifying 
a point in this space specifies a state of the system, and vice 
versa. However, in order to apply the concept of phase space to 
a time series, which is a scalar sequence of measurements, one 
has to convert it into a set of state vectors. 

For example, let us suppose we have an observed time se¬ 
ries, x{ti), with i = At,2At, ...,NAt. One can construct a set 
of 2(3, 4, ...)-dimensional vectors by selecting pairs (triplets, 
quadruplets, ...) of data points, whose count rate values rep¬ 
resent the values of the different components of each vector. 
The data points defining the vector components are chosen such 
that the second point is separated from the first by a time de¬ 
lay At', while the third is separated from the first by 2At', and 
so on. This process, called time-delay reconstruction, can be 
easily generalized to any n-dimensional vector. The time de¬ 
lay At' is usually much larger than the sampling interval At of 
the observed light curve. In fact, in non-linear dynamics, where 
this technique is used to determine fractal dimensions and dis¬ 
criminate chaos from stochasticity, the choice of the time delay 
is not arbitrary, but should satisfy specific rules (see Kantz & 
Schreiber 1997 for a review). 

The set of n-dimensional vectors, derived from the orig¬ 
inal light curve, is then mapped (or embedded) into an n- 
dimensional vector space, producing a phase space portrait, 
whose topological properties represent the temporal properties 
of the corresponding time series. 

In order to quantify the characteristics of a phase space 
portrait, a widely used method is the correlation integral C{R) 
(Grassberger & Procaccia 1983). An exhaustive description of 
the correlation integral method and its application to X-ray 
light curves of AGN is given by Lehto et al. (1993). 

Another method that can quantify the topological proper¬ 
ties of a phase space portrait is the scaling index method (here¬ 
after SIM; e.g., Atmanspacher et al. 1989) which has been em¬ 
ployed in a number of different fields (e.g.. Rath & Morfill 
1997; Rath et al. 2002), including the study of AGN time series 
in the case of Ark 564 (Gliozzi et al. 2002). Differently from 
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Fig. 7. Scaling index distributions of 3C 390.3 during the 1999 
(solid line) and 2000 (dashed line) monitoring campaign, us¬ 
ing light curves binned at 32 s (top panel), and 160 s (bottom 
panel). The thick continuous lines represent the scaling index 
mean value during 1999, {Q' 99 ), while the thick dashed lines 
represent (ago)- 

the correlation integral, which is a global indicator, the scaling 
index method measures both the global and local properties of 
a phase space portrait. Indeed SIM is based on the local es¬ 
timate of the correlation integral for each vector in the phase 
space, by measuring the “crowding” of the vectors around each 
individual point in the phase space. 

For example, suppose we choose to define vectors of di¬ 
mension D from a time series x(f,), with time delay At' = At. If 
the time series has N observed points, the number of vectors is 
N' = N - D + 1. Then, for each of the N' vectors, say the /-th, 
we calculate the cumulative number function, 

Q(/?) = n{Mj < n (1) 

In other words, Ci(/^) measures the number of vectors J 
whose distance dij from the vector i is smaller than /?. For es¬ 
timating the local scaling properties of the distribution of vec¬ 
tors, one approximates the function C,(f?) with a power law of 
the form C,(/^) ~ within a given interval [/?i,/? 2 ] (which is 
related to the typical distances between the data points, that, in 
turn, depend on the choice of the embedding space dimension). 
The exponent a, is called scaling index.ln practice, the scaling 


index for each vector i is computed by considering the logarith¬ 
mic increase of Ci{R) within a suitably chosen interval [f?i, /? 2 ]- 
In all our calculations a, has been determined as follows: 

logC,(/?i,/? 2 )-logC,(/?i) 

O'; = -■ (2) 

(log /?2 -log/?l) 

where C,(f?i,f? 2 ) denotes the mean cumulative number within 
the interval [f?i,/? 2 ]- Note that, by definition, the scaling index 
cannot have negative values. 

The scaling index can be estimated for all the N' vectors 
in the phase space, producing a distribution of scaling indices, 
whose shape, width, and centroid quantify the the topological 
properties of the phase space portrait. In this way, the tempo¬ 
ral properties of the original light curve are translated into the 
properties of the scaling index distribution. 

Note that for a purely random process (e.g., a process with 
an associated time series where the points vary independently) 
the shape of the histogram is broad, dominated by the right- 
ward tail, and the average scaling index {a) tends to the value 
of the dimension of the embedding space. On the other hand, 
for highly correlated processes and for deterministic (chaotic) 
processes, the histogram is typically narrowly shaped, with the 
value of {a) always significantly smaller than the dimension of 
the embedding space. 

One of the most appealing characteristics of SIM is its abil¬ 
ity to discern underlying structures in noisy data. Therefore, it 
might be more sensitive than linear methods like PSD, SF, and 
PDF, in identifying intrinsic differences in the temporal prop¬ 
erties of 3C 390.3 during 1999 and 2000. 

To test the possible presence of two distinct temporal states 
of the system in the 1999 and 2000 observations, we have es¬ 
timated the distribution of scaling indices of these two light 
curves. We stress that, at this point, we are not interested in 
determining any kind of “dimension” of the system, since this 
would be meaningless for non-deterministic systems. Our aim 
is to compare the 1999 and 2000 histograms of the scaling in¬ 
dex values, and try to detect any significant differences between 
them. To the extend that these histograms correspond to the 
temporal properties of the system, any differences will imply 
variations in the intrinsic temporal properties of the system. 

We calculated the scaling index for all the vectors of the 
two phase space portraits derived from the two monitoring 
campaigns, using embedding spaces of dimensions ranging be¬ 
tween D — 2 and 10. There are no systematic prescriptions 
for the choice of the embedding dimension. The discriminating 
power of the statistic based on the scaling index is enhanced 
using high embedding dimensions. This can be understood in 
the following way: if the data are embedded in a low-dimension 
phase-space many of them will fall on the same position, los¬ 
ing part of the information. On the other hand, the choice of 
too high embedding dimensions would reduce the point density 
in the pseudo phase-space, lowering the statistical significance 
of the test. As for the time delay At', we have used At' =32 
s and 160 s, producing, with the former choice, ~4500-5000 
vectors, and, with At' =160 s, -1000-1200 vectors, depending 
on the dimension and on the campaign (the 2000 campaign, 
albeit shorter, has longer individual observations, producing a 
slightly higher number of data points and thus of vectors). It 
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Fig. 8. Top left panel: Distribution of Aasynth obtained from red-noise simulations (tbin= 160 s) compared with Aoobs (thick solid 
line). The dotted, dashed and dot-dashed lines represent the 95%, 99%, and 99.9% confidence levels, respectively. Top middle 
and right panels: Distribution of mean synthetic scaling indices obtained from red-noise simulations, (aggiooisynth), compared 
with (Q; 99 (oo),ois) (thick solid line). Bottom left panel: Distribution of AcTsynthi i-C-, of the differences between standard deviations 
around (a) of the synthetic distributions of scaling indices. Bottom middle and right panels: Distribution of synthetic standard 
deviations around (a) in 1999 and 2000, respectively. 


is worth noticing that the use of small time bins cause the data 
points to be correlated, but this is not necessarily a drawback, 
in particular in the case where one wants to measure a differ¬ 
ence in the degree of correlation of two time-series. In all the 
cases (i.e. for all D values) we found that the 1999 scaling in¬ 
dex histogram is systematically different from the 2000 his¬ 
togram, with the former being significantly narrower, and its 
centroid being smaller than the latter. This result suggests that 
the points in the 1999 light curve are more “strongly” corre¬ 
lated (e.g., over longer time-scales) than those in the 2000 light 
curve. 

Importantly, the results obtained using tbin=32 s are fully 
consistent with those obtained with tbin=160 s, suggesting that 
the scaling index method is actually able to detect the same un¬ 
derlying structure despite the different levels of Poisson noise 
(by using tbin=32 s, the ratio between the mean error and 
the mean count rate {err) j{rate) increases by a factor larger 
than 80% compared to tbin=160 s). The resulting histograms 
for a 5-dimensional phase space for tbin=32 s (top panel) and 
tbin=160s (bottom panel) are plotted in Fig.^l the histograms 
with solid lines represent data from the first monitoring cam¬ 


paign, whereas the histograms with dashed lines refer to 2000. 
In both panels, the thick solid (dashed) line represents the lo¬ 
cation of the centroid of (a) in 1999 (2000), which remains 
roughly at the same location despite the large change in the 
signal-to-noise level. 

4.4.1. Impact of red and Poisson-noise effects on the 
SIM 

In the previous sections we have have seen how the red-noise 
effects might affect the results from the SF and PDF analyses, 
casting doubts on the statistical significance of apparently re¬ 
markable differences observed in the temporal properties be¬ 
tween 1999 and 2000. Therefore, it is imperative to assess the 
impact of the red-noise effects on SIM. 

Before discussing the results of red-noise simulations, 
it is important to note that in May 2000 (i.e., near the be¬ 
ginning of the second monitoring campaign of 3C 390.3) 
the propane layer on PCUO was damaged, causing a sys¬ 
tematic increase of the background. This effect is accounted 
by the most recent PCA background models and does not 
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Fig. 9. Left panel; Distribution of Aasymh obtained from red-noise simulations with 50% extra noise compared with Aoobs indi¬ 
cated by the thick continuous line (tbin= 160 s). Middle and right panels; Distribution of mean synthetic scaling indices obtained 
from red-noise simulations with 50% extra noise, {agg(oo)^synth, compared with {a' 99 (oo),oij (thick continuous line). The dotted and 
dashed lines represent the 95% and 99% confidence levels, respectively. 


represent a problem when sufficiently large time bins are used 
(http;//lheawww.gsfc.nasa.gov/users/craigm/pca-bkg/bkg-users 
However, the SIM, being a statistical method, needs large 
amounts of points to work properly and thus necessarily re¬ 
quires the use of small time bins. We therefore paid particular 
attention at possible systematic effects caused by the use of 
PCUO data in 2000. Indeed, the SIM does notice a difference 
between the PCUO and PCU2 during 2000, with the former 
PCU being characterized by larger values of (a). To be on the 
safe side, we report the SIM results obtained excluding the data 
from PCUO in 2000. However, for the sake of completeness, 
we also carried out all tests using both PCUs in 2000, and 
reached the same conclusions. 

To assess the impact of red-noise effects on SIM, we have 
used 10^ synthetic light curves obtained from Monte Carlo 
simulations very similar to those described in §4.2, the only 
difference being the use of one PCU (PCU2) instead of two 
(PCU0-I-PCU2) for the simulated data of 2000. As before, we 
report our results in the case when the PSD break at 1/300 
day“*. 

For each pair of synthetic light curves the SIM is applied 
and the difference Aa = (aoo) - {agg) is computed. For sim¬ 
plicity and computational reasons (since we are dealing with 


thousands of simulations of data-sets containing thousands of 
htfMlInts), we limit our discussion mainly to Aa and Actq, (i.e., 
the difference between the dispersion of the a values around 
(a)) as the discriminating parameters between the scaling in¬ 
dex histograms of the two light curves. However, it must be 
kept in mind that SIM provides more information, and one can 
use the full distribution of the scaling index values for the most 
accurate comparison between the two observed histograms. 

Fig. 13 (left panels) reveals that red-noise effects are unable 
to explain Ao'ob.s (top panel; thick solid line) nor AcTobs (bottom 
panel) at very high significance level (P< 0.1%). It is worth 
noting that, if no Poisson noise was added to the synthetic light 
curves, then the distributions of Aa and Act would be centered 
around zero. The fact that the distributions are shifted towards 
positive values is due to the use of two PCUs in 1999 and only 
one in 2000, and hence to the different Poisson noise contribu¬ 
tions to the 1999 and 2000 synthetic data. 

To understand better these results, we have also plotted the 
distributions of the scaling index centroids and standard devia¬ 
tions in 1999 and 2000 (Fig.|3middle and right panels, respec¬ 
tively). It is clearly shown that the distributions of simulated 
data of 1999 (Fig. |3 middle panels) are marginally consistent 
with the corresponding observed values; the chance probabili- 
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ties to obtain values consistent with (agg) and cr^^gg are P~ 6% 
and P~ 15%, respectively. In contrast, the distributions of sim¬ 
ulated data of 2000 (Fig. [fright panels) are totally inconsistent 
with the observed values (P< 0.1% for both (aoo) and cr^^o 
distributions). 

The fact that AcTsynth and Acrsy„th are not centered around 
zero, because of the presence of different Poisson noise com¬ 
ponents in the two light curves, motivated us to investigate fur¬ 
ther the the impact of the Poisson noise on the distribution of 
the scaling index values. In principle, the effects of Poisson 
noise have been already accounted by the red-noise simula¬ 
tions, where the “appropriate” noise level has been computed 
also considering the different background levels during 1999 
and 2000. Nevertheless, we considered what the effects on the 
SIM results would be if we add, arbitrarily, some artificial ex¬ 
tra Poisson noise (equal to half of the noise actually present 
in the data) to the synthetic light curves. In other words, we 
consider the possibility that the “real” instrumental noise in the 
light curves is, for some unknown reason, underestimated in 
the previous simulations, and what effects this may have on the 
SIM results. 

The results of the noisy simulations are plotted in Fig. |9l 
Based on the Aasynth distribution only (top left panel), we 
would conclude that the simulated results are in agreement 
with the observed value of Aa (the chance probability is indeed 
P~ 17%). This might lead to the erroneous conclusion that the 
differences observed using SIM could be simply ascribed to the 
fact that the Poisson noise in the previous simulations was not 
accounted for correctly. However, a look at the AcTsynth distri¬ 
bution (Fig. 1^ bottom left panel; P< 0.1%) reveals that this 
is not the case. The same conclusion is reached looking at the 
synthetic distributions in 1999 and 2000, and keeping in mind 
that, in order to infer consistency between simulated and ob¬ 
served values, the distributions of both scaling index centroids 
and standard deviations should be consistent with the respec¬ 
tive observed values. This is clearly not the case in 1999 and in 
2000 ; in the first case, the synthetic distribution of {agg) is very 
poorly consistent with the observed value (P= 0.1%), whereas 
in 2000, it is the synthetic distribution of 0 -^^ values which is 
entirely inconsistent with (P< 0.1%). 

Another macroscopic consequence of the arbitrary increase 
of the Poisson noise is the the substantial right-ward shift of 
(a). This might appear at odds with the lack of shift in (a) ob¬ 
served in the real data when passing from tbin= 160 s to tbin=32 
s (see Fig.0. However, this apparent discrepancy can be eas¬ 
ily explained by the most important property of the SIM, i.e. 
its ability to discern an underlying structure in noisy data, if an 
underlying structure is present. On the other hand, if no under¬ 
lying structure is present (as in the case of the simulated red- 
noise data), the increase of noise (and hence of randomness) 
simply causes a shift of (a) to larger values. 

An alternative and model-independent way to assess if 
Affiobs is dominated by Poisson noise effects is based on the 
shuffling of the temporal order of the RXTE orbits. In this 
way, the original long-term structure is destroyed, while the 
Poisson noise levels and all the statistical moments of the dis¬ 
tributions in 1999 and 2000 are unaffected. If the orbit-shuffled 
light curves yield a Aa consistent with Aaobs, then we should 


conclude that the difference in Poisson noise between the 1999 
and 2000 (along with the difference in the short-term intra¬ 
orbit variability) is the simplest explanation of the difference 
in scaling index observed. On the other hand, if the values of 
Aashuffled are not consistent with the observed value, the direct 
consequence is that the long-term temporal structure of the two 
original light curves does play an important role. 

The results of this test, carried out using light curves with 
tbin=32 s because of their higher Poisson noise level (but sim¬ 
ilar results are obtained with tbin=160 s), are shown in Fig. Ci 
The left panel clearly shows that the orbit-shuffled results are 
inconsistent with the Aaob.s at very high significance level. The 
middle and right panels show the distributions of the centroid 
(a) and the width cr„ of the scaling index histograms obtained 
by shuffling the orbits order in 1999 (very similar results are ob¬ 
tained with 2000 data). Both plots show that the orbit-shuffled 
quantities narrowly cluster at a location, which is, statistically 
speaking, far away from the real data, indicating that the long¬ 
term temporal structure contributes substantially to the SIM re¬ 
sults. 

We conclude that, contrary to the results based on the com¬ 
parison between the 1999 and 2000 SFs and PDFs, the ob¬ 
served differences between the respective distributions of the 
scaling index values are almost certainly intrinsic, i.e. they can¬ 
not be explained by Poisson noise, red-noise, nor any other in¬ 
strumental effect that we can think of. To the extend that the 
differences of the scaling index histograms correspond to dif¬ 
ferences in the temporal properties of the two light curves, we 
consider them as significant evidence of non-stationarity in the 
the X-ray emitting process in 3C 390.3. 

5. Spectral Variability 

After searching systematically for the evidence of temporal 
changes in 3C 390.3, we now want to investigate whether the 
apparent changes of the statistical timing properties are accom¬ 
panied by spectral variations, and what is the nature of the spec¬ 
tral changes (e.g., is 3C 390.3 experiencing an actual state tran¬ 
sition?). To this end, we have performed two different kinds of 
analysis. 

1) The first (model-independent) method is based on the X- 
ray color - count rate plot, often used for GBHs to separate the 
spectral states. The plot of X-ray color (i.e., the ratio of the 7- 
20 over the 2-5 keV count rate) versus the total 2-20 keV count 
rate is shown in Fig. CH The hardness ratio clearly decreases 
with increasing flux, during the second part of the 1999 and the 
whole of the 2000 campaigns. This trend is typical of type 1 
radio-quiet AGN, and has been observed many times in the past 
(e.g., Papadakis et al. 2002; Taylor et al. 2004, and references 
therein). This supports the conclusions of Gliozzi et al. (2002), 
that the X-ray emission in 3C 390.3 is not jet-dominated. 

During the first half of the 1999 campaign, the source was 
at its lowest flux state. The spectrum of the source also appears 
to be quite hard, consistent with the general trend of spectral 
hardening towards low flux states. However, the scatter of the 
hardness ratio points is very large, and it seems that during this 
state, the spectral variations are more or less independent of the 
source flux variations. So, based on the color-count rate plot. 
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Fig. 10. Left panel; Distribution of Acr obtained by shuffling the time order of the RXTE orbits. The thick solid line is the location 
of Aaobs- Middle panel; Distribution of (a) (short, solid lines) obtained by shuffling the orbits order in 1999 (tbin=32 s). The 
open diamond gives the mean value of the shuffled data and its error bar represents the 3cr confidence level. The long, dashed line 
represents the original value, obtained without any shuffling. The same symbol convention is used for the right panel, showing 
the distribution of ctq. in 1999. 
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Fig. 11. 7 -20 keV/2-5 keV X-ray color plotted versus the total 
count rate of 3C 390.3. Star symbols correspond to the hrst 
200 days of the 1999 light curve, the hlled dots describe the 
second part of the 1999 campaign, and the open diamond the 
2000 light curve. For the sake of clarity, only a characteristic 
error for each cluster of points has been plotted. All the errors 
associated with the X-ray colors are reported in Fig.^bottom 
panel. Time bins are 5760s (~ 1 RXTE orbit). 


it seems possible that during this period the X-ray emission 
mechanism may behave in a different way, although the errors 
associated with the X-ray colors (see Fig.^bottom panel) pre¬ 
vent us from reaching a hrm conclusion. 

2) More stringent constraints can be derived by using a second 
(model-dependent) method; the time-resolved spectral analy¬ 
sis. Since the data consist of short snapshots spanning a long 
temporal baseline, they are well suited for monitoring the spec¬ 
tral variability of the source. 

We divided the light curve of Figure ^ into thirteen inter¬ 
vals, of duration 50-60 days each, with net exposures times 
ranging between 15 ks and 32 ks. This scheme represents a 
trade-off between the necessity to isolate intervals with differ- 


Table 1. Time Intervals Used for Time-Resolved Spectral 
Analysis 


Start Time 

End Time 

Exposure 

Flux" 

r 

(y/m/d h;m) 

(y/m/d h;m) 

(ks) 



99/01/08 00;26 

99/02/25 17;58 

17.90 

1.87 

1 55^-^^ 

'^'^ 0.05 

99/02/28 16;48 

99/05/11 01;05 

24.10 

1.68 

1 54OO3 

^•*^^ 0.05 

99/05/14 00;58 

99/07/10 20; 15 

21.20 

2.64 

i-KS 

99/07/13 08;56 

99/09/11 16;03 

22.56 

3.86 

1 - 6 C 3 

99/09/14 18;07 

99/11/13 22;58 

23.38 

3.71 

^•^'0.02 

99/11/17 14;25 

00/01/03 08;00 

18.69 

4.71 

1 740.02 
^ •' ^ 0.02 

00/01/06 11;50 

00/02/29 06;57 

15.26 

3.84 

^•^^ 0.03 

00/03/03 04;26 

00/05/11 12;43 

32.46 

4.90 

^ •' ^ 0.02 

00/05/14 08;41 

00/07/04 01;22 

22.50 

3.80 

^•^^0.02 

00/07/07 00;09 

00/09/01 22;29 

26.91 

4.10 

^•^‘^0.02 

00/09/05 5;52 

00 / 11/10 12 ; 12 

32.22 

3.91 

^‘^^0.02 

00/11/13 12;57 

00/12/31 03;12 

25.47 

5.03 

1 67®-^^ 
^•^' 0.02 

01/01/03 08;53 

01/02/24 00;26 

24.45 

3.46 

^•^^0.02 


“ 2-10 keV flux in units of (10 "erg cm ^ s *)> calculated assuming a 
power law plus Gaussian line model. 


ent but well dehned average count rates and the need to reach a 
signal-to-noise ratio (S/N) high enough for a meaningful spec¬ 
tral analysis. The dates, exposure times and mean fluxes in the 
thirteen selected intervals are listed in Table ^ 

We fitted the 4-20 keV spectra with a model consisting of 
a power-law, absorbed by the Galactic absorbing column, plus 
a Gaussian Fe Ka line with a rest energy fixed at 6.4 keV. In all 
the intervals, the line is statistically required at more than 99% 
confidence level, according to an F-test. However, due to the 
low S/N, the line parameters are characterized by large errors 
(~ 100% on the line flux and width) and thus their evolution 
in time cannot be investigated properly. On the other hand, the 
photon indices are well constrained with uncertainties of the 
order of 1-3%. Figure El shows the plot of the photon index 
F versus the 2-10 keV flux, with the filled dots representing 
intervals during the 1999 campaign and open diamonds during 
the 2000 campaign. The photon index increases steadily with 
the 2-10 keV flux, without any indication of a discontinuity be¬ 
tween the F corresponding to the low state and the other photon 
indices. 
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fbreak moviiig from ~ 1/300 - 1/400 days to ~ 1/50 -1/100 days. 
However, this apparent shift is not of high statistical signifi¬ 
cance. The Structure Function analysis also indicates a similar 
change in the characteristic time-scale of the system. Both SFs 
have a power-law shape up to lags shorter than ~ 60 days. At 
longer time-scales, the 2000 SF appears to flatten to a roughly 
constant level, while the 1999 SF keeps rising up to the largest 
time lags sampled. Finally, the comparison between the two 
Probability Density Functions (see Fig. |3, suggest the exis¬ 
tence of an intrinsic difference between the two time series, 
with the 1999 data characterized by a bimodal distribution and 
the 2000 data following a uniform distribution. 

However, based on the results from Monte Carlo red-noise 
simulations, which were employed in order to test the signifi¬ 
cance of the apparent differences between the 1999 and 2000 
SF and PDFs, we found that they can be the result of the red- 
noise character of the variability process with a probability at 
least as high as ~ 10%. 


Fig. 12. Photon index, F, plotted versus the 2-10 keV flux 
(10 ergcm^^ s ')- Filled dots refers to the 1999 campaign, 
and the open diamond to the 2000 light curve. The dashed line 
represents a linear least squares fit to all the data points. 

Based on this analysis, we conclude that 3C 390.3 is not un¬ 
dergoing a transition between two distinct spectral states during 
the RXTE ’99 and ’00 monitoring campaigns. 

6. Discussion 

We have studied two one-year long, well sampled X-ray light 
curves of 3C 390.3 from the RXTE archive, in search of ev¬ 
idence for variations in its temporal properties. Our motiva¬ 
tion stems from the fact that a long-term ROSAT light curve of 
the this broad-line radio galaxy had shown evidence for non¬ 
linearity and possibly non-stationarity in the past (Leighly & 
O’Brien, 1997). The light curves we have used are the longest, 
and best sampled light curves on time-scales of days/months 
of this source to date. Hence, they offer a good opportunity to 
investigate the important issue of stationarity in the X-ray light 
curves of AON. 

6.1. Multi-technique Timing Approach: Linear Methods 

In order to assess whether the light curve of 3C 390.3 does 
show variations in its temporal properties, we performed a thor¬ 
ough temporal analysis based on different complementary tim¬ 
ing techniques and extensive Monte Carlo simulations. The 
first advantage of this approach is that it leads to more ro¬ 
bust results, not biased by the specific characteristics of a sin¬ 
gle method. In addition, since the diverse methods employed 
probe different statistical properties, they provide us with com¬ 
plementary pieces of information that can be combined to con¬ 
strain the physical mechanism underlying the observed vari¬ 
ability. 

The Power Spectral Density analysis results suggest a pos¬ 
sible shift of a frequency break between 1999 and 2000, with 


6.2. The Scaling Index Method 

In contrast with the linear methods, that showed only sugges¬ 
tive evidence of a change in the temporal properties of 3C 390.3 
but were unable to rule out a major influence of red-noise ef¬ 
fects, the Scaling Index method clearly reveals the presence 
of intrinsic changes that cannot by any means be explained by 
red-noise nor Poisson effects (P<0.1%). 

It is difficult to understand the consequences of this result 
at this point. One obvious possibility is that the variability pro¬ 
cess of 3C 390.3 holds more information than the red-noise 
PSD reveals. The main property of the synthetic light curves 
that we have produced using the method of Timmer & Koening 
(1995) is that their PSD has a specific red-noise shape. The 
disagreement then between the results from the simulated light 
curves and the observed SIM values could imply that there is 
“something” more than just “red plus Poisson noise” in the ob¬ 
served light curves. In fact, there have been indications the last 
few years that this may well be the case. For example, Uttley, 
McHardy & Vaughan (2004) have shown that the rms-flux re¬ 
lation recently discovered in the X-ray light curves of AGN 
and GBHs implies that their light curves are intrinsically “non¬ 
linear”. The simulations that we have used in this work to test 
the SIM results do not reproduce this rms-flux relation. Perhaps 
then, the X-ray light curve of 3C 390.3 is non-linear as well, 
and the differences that we detect with the SIM method re¬ 
flect this property. Advanced simulations (which are beyond 
the scope of the present work) are needed in order to inves¬ 
tigate this issue further, and examine whether non-linear syn¬ 
thetic light curves which do reproduce the rms-flux relation are 
consistent with our results. 

On the other hand, it is also possible that, irrespective of 
whether the light curves are linear or non-linear, the SIM dif¬ 
ferences that we observe imply real differences in the temporal 
properties of the source between 1999 and 2000. The impor¬ 
tant issue of course is to identify these properties, in order to 
advance our understanding of the variability mechanism. This 
is not an easy task, because the nature of the AGN time variabil- 
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ity is still poorly known (apart from the “linear vs. non-linear” 
issue, one has to consider other issues as well like “stochas¬ 
tic vs. chaotic” for example). At the moment, no clear physical 
insights can be derived from this non-linear tool. Indeed, as 
pointed out by Vio et al. (1992), a naive application and inter¬ 
pretation of the non-linear methods results may lead to concep¬ 
tual errors. 

Nevertheless, the much narrower 1999 scaling index dis¬ 
tribution, with its centroid (a) located at a value significantly 
lower than (aoo), is consistent with the assumption that the 
1999 time series is characterized by a “higher degree of correla¬ 
tion”. Since the correlation structure of a time series is charac¬ 
terized by its auto-covariance function, perhaps this result does 
suggest that the the width of the auto-covariance function has 
decreased during the 2000 observations. This is qualitatively 
in agreement, with the hypothesis that the characteristic time- 
scale during 1999 is longer than during 2000, as suggested by 
PSD and SF analyses. However, at this point this is merely a 
suggestion and alternative hypotheses cannot be ruled out. If 
true, the differences in the scaling indices do not challenge the 
“red plus Poisson noise” view; they simply demonstrate that 
the first term in this expression may vary with time. 

6.3. 3C 390.3 and the analogy between AGN and 
GBHs 

It has been suggested many times in the past, and is now widely 
believed, that AGN are simply the high mass analogs of GBHs. 
For example, PSD analyses of long, high signal-to-noise, well 
sampled X-ray light curves of a few AGN have shown clearly 
that their power spectra are very similar to those of GBHs (e.g. 
Uttley et al., 2002; Markowitz et al., 2003; McHardy et al., 
2004). 

One of the most interesting aspects of the X-ray variabil¬ 
ity behavior of the GBHs is that they often show transitions to 
different “states”. “State” transitions have not been observed so 
far in AGN. The energy spectral results reported in Section 5 
clearly argue against the presence of a bona fide state transi¬ 
tion in 3C 390.3 as well. This is not unexpected if we consider 
the time-scales at play: in GBHs state transitions take place on 
time-scales of days (e.g. Zdziarski et al. 2004) which translate 
into thousands of years for 3C 390.3, assuming a linear scaling 
between the black hole masses. 

However, even within the same state, the timing proper¬ 
ties of GBHs, like the PSD characteristic time-scales, vary with 
time. For example, during the “low/hard state”, the characteris¬ 
tic time-scales in the X-ray power spectrum of Cyg X-1 change 
with time in a correlated way with the spectral variations of the 
source (Belloni & Hasinger, 1990; Pottschmidt et al., 2003). 
This is one example of “stationarity-loss”, which can offer im¬ 
portant clues on the mechanism responsible for the observed 
variations. These variations of the temporal properties, can take 
place on time-scales as short as tens of seconds (e.g. Casella et 
al. 2004). These time-scales are consistent with the length of 
the 1999 and 2000 observations of 3C 390.3 which we study in 
this work, assuming again, that characteristic time-scales scale 
linearly with black hole mass. 


We find that the energy spectrum of the source in the first 
part of the 1999 light curve is significantly harder than that of 
the high flux peaks after day ~ 400 (Fig. 1). Furthermore, the 
SIM results shows variations in the timing properties of the 
source during the same periods. In itself, this result strengthens 
the analogy between the X-ray variability properties of AGN 
and GBHs: just like in the Galactic systems, 3C 390.3 does 
show simultaneous spectral and intrinsic timing property vari¬ 
ations. 

However, as we discussed above, the SIM results are not 
so easy to interpret at the moment. On the other hand, the re¬ 
sults presented in sections 4.1-4.3, do suggest a much closer 
analogy between 3C 390.3 and say Cyg X-1, but they are not 
statistically significant, possibly because the quality of the light 
curves at hand is not good enough. 

Our results show that the spectral slope, F, was ~ 1.55 dur¬ 
ing the first ~ 150 days of the 1999 monitoring observations, 
increasing by a factor of ~ 10%, during the peak flux periods 
of the 2000 monitoring observations. The Cyg X-1 slope does 
not become as hard as 1.55, nevertheless, according to the re¬ 
sults of Pottschmidt et al. (2003), when F increases from its 
lowest value of 1.8 to a value of ~ 1.95-2 (an increase by a 
factor of ~ 10% as well) the centroid frequency of the L 2 com¬ 
ponent in its power spectrum increases by a factor of ~ 5 - 6 
(as shown in the middle panel of their Fig. 7). Since the L 2 
centroid frequency is commonly thought to correspond to the 
“high-frequency” PSD breaks in AGN, this is the same behav¬ 
ior that the results of 3C 390.3 suggest. 

As mentioned above, we believe that the temporal and spec¬ 
tral behavior of 3C 390.3, as seen by RXTE during the 1999 and 
2000 monitoring campaigns, qualitatively mimics the proper¬ 
ties of the GBH Cyg X-1 during the low/hard state, provided 
that the appropriate time-scale differences are taken into ac¬ 
count. 

7. Summary and Conclusions 

The main results and conclusions of this work can be summa¬ 
rized as follows: 

- We have studied in detail the variability characteristics of 
the two-year long, RXTE monitoring light curve of 3C 
390.3. Linear methods, like PSD, SF and PDF analysis 
show suggestive evidence of a change in the characteristic 
time-scale of the system between 1999 and 2000. 

- However, only the results from the scaling index method 
show a clear and significant evidence of a change in the 
intrinsic properties of the system that cannot be explained 
by red or Poisson noise effects. 

- It is not easy to interpret the SIM results at this moment. 
To the extend that the scaling index is related with the tim¬ 
ing properties of the variability process, one possibility is a 
change of the characteristic time-scale between 1999 and 
2000 , in agreement with the suggestive results from the 
PSD, SF and PDF analyses. 

- Our results reinforce the similarity between the X-ray vari¬ 
ability properties of AGN and GBHs. In analogy with GBH 
phenomenology, the loss of stationarity in 3C 390.3 seems 
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to correspond to an increase in the frequency break as the 
spectrum softens, as observed in Cyg X-1, when the system 
is in its low/hard spectral state. 

- Since the SIM is not yet a standard technique, further work 
is needed in order to investigate its properties and its possi¬ 
ble relation with more frequently used (and better under¬ 
stood) statistical analysis tools like the power spectrum. 
Nevertheless, we believe that the present results demon¬ 
strate the benefits from the application of the SIM (and pos¬ 
sibly other non-linear methods) in the analysis of the AGN 
X-ray light curves. 

- The RXTE archive is populated with numerous, long AGN 
light curves. This is a good opportunity to start a system¬ 
atic study of these data, for a search of spectrally related 
variations of the timing properties of these systems. Since 
the outcome of any Monte Carlo simulation based on an 
assumed “model-PSD” is inevitably dependent on the as¬ 
sumptions made, we believe that a thorough analysis of a 
large number of high-quality AGN light curves is the most 
direct way to examine whether the present results are in¬ 
dicative of the AGN behavior in general, or just a peculiar¬ 
ity associated with 3C 390.3. We plan to explore this issue 
in a future work. 
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